Anomalies in the NMR of Silicon: Unexpected Spin Echoes in a Dilute Dipolar Solid 
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NMR spin echo measurements of ^^Si in Silicon powders have uncovered a variety of surprising 
phenomena that appear to be independent of doping. These surprises include long tails and even-odd 
asymmetry in Carr-Purcell-Meiboom-Gill (CPMG) echo trains, and anomalous stimulated echoes 
with several peculiar characteristics. Given the simplicity of this spin system, these results, which 
to date defy explanation, present a new and interesting puzzle in solid state NMR. 



cn 
o 
o 

< 



> 

oo 
in 

00 

o 
m 
o 



o 



PACS numbers: 03.65.Yz, 03.67.Lx, 76.20.+q, 76.60.Lz 



In order to implement quantum computation (QC) 
based upon spins in semiconductors 0,13, 0, S Hi ^ 
detailed understanding of spin dynamics in these mate- 
rials is required. To this end, we carried out a series 
of NMR measurements that were motivated by a sim- 
ple question: what is the ^^Si decoherence time {T2) in 
Silicon? Earlier NMR studies in Silicon addressed other 
questions |EI3,I3- 

We find that it is possible to detect the ^ssi (4.67% 
natural abundance (n.a.), spin-i) NMR signals out to 
much longer times than was previously thought possible, 
and so far, we have been unable to explain these results in 
terms of well-known NMR theory Surprises in such 
a simple spin system appear brand new to NMR, and 
understanding their origin is of fundamental importance. 
In this paper, we describe the phenomena and recount 
tests we have made to explore possible explanations. 

Two standard experiments that measure T2 are re- 
ported. First, using the Hahn echo sequence (HE: 90x- 
(^)-180y-(^)-ECHO 113), the measured decay, with 
T2he ~ 5.6 msec, is in quantitative agreement with 
that expected for the static ^^Si-^^Si dipolar interac- 
tion. This decay mechanism is commonly encountered 
in solids, and a number of ingenious pulse sequences 
have been invented to manipulate the interaction Hamil- 
tonianjpushing echoes out to times well beyond T2„j^ 
|lIillilElIllElEIi3-A common thread running 
through those sequences is the use of multiple 90° pulses, 
and pulses applied frequently compared to r2^£, , which 
refocus the homonuclear dipolar coupling. The same 
cannot be said about the second sequence that we used 
to measure T2, the Carr-Purcell-Meiboom-Gill sequence 
(CPMG: 90x-{{^) -180f-(^)-ECHO}''^p'='^* 
[13 )• Specifically, the CPMG sequence is not expected to 
excite echoes beyond T2hj^ , since 180° pulses should not 
affect the bilinear homonuclear interaction. This state- 
ment is exact in two important limits: either for unlike 
spins or for magnetically-equivalent spins. 

Therefore, we were surprised to find that CPMG echoes 
are detectable long after T2he) '^^'^ echo peaks ap- 
pear nearly identical in Silicon samples with very dif- 
ferent dopings. This CPMG "tail" appears to be even 
larger at low temperatures. In addition, as the interpulse 
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a=Si:B N.=10 cm 0.37 kHz 



8 T, 
0.66 cm 2.8 hr 



b=Si:Sb No=10 cm 0.32 kHz 0.14 cm 4.6 min 
c=Si:P ND=10^^cm"^ 0.35 kHz 6.40 cm 5.5 hr 
3.17 kHz 0.03 cm 4.8 sec 
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FIG. 1: Two standard measurements of the ^^Si T2 in pow- 
dered Silicon at room temperature (RT). CPMG echo trains 
are shown for four samples (a-d) with different doping, full 
width half max (FWHM) , skin depth (S) , and spin- lattice re- 
laxation time (Ti). Since samples a, b, and c exhibit much 
wider echoes than sample d, only the top portion of their 
echoes are visible. Hahn echo measurements (circles with 
crosses for sample d, others are supressed for clarity) agree 
quantitatively with the dipolar decay curve (solid line) cal- 
culated for the Silicon lattice (Eqn. 4, see text). Despite 
big changes in doping (e.g., xlO® in P-concentration between 
samples c and d), the peaks of the CPMG echoes are nearly 
identical to each other, and they are detectable long after 
the Hahn echoes decay to zero. These measurements are in 
a 7.027 Tesla field {B \\ z, with = 59.48 MHz). N.B., ah 
figures are in color online. 



spacing (TE) is increased, the CPMG echoes develop a 
pronounced "even-odd asymmetry" (e.g., long after spin 
echo #1 ("SEl") is in the noise, spin echo #2 ("SE2") is 
clearly observable). Lastly, we show how an "anomalous 
stimulated echo" is observed in this system, with several 
peculiar characteristics. 

Figure 1 shows CPMG echo trains acquired in four 
different Silicon samples (both n-type and p-type). 
As the legend shows, the ^^Si NMR spectrum (0.3 
kHz<FWHM<3 kHz), the echo shape, and the Ti (from 
4.8 sec to 5.5 hours at RT) can be quite different, for sam- 
ples with wide variations in doping 0|. Despite these 
big changes (e.g., xlO^ in P-concentration), the peaks of 
the CPMG echoes are nearly identical to each other, and 
they persist long after the Hahn echoes have died away. 
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Qualitatively, the long tail evokes a well-known effect 
in liquid-state NMR, where diffusion causes slow changes 
in the local field leading to an extrinsic decay of the Hahn 
echoes 0. Applying frequent refocusing pulses renders 
the dynamics "quasi-static" , enabling the CPMG echoes 
to persist to longer times, and revealing the intrinsic T2. 
However, in our data, the Hahn echoes appear to persist 
out to the "intrinsic" T2 curve, and the CPMG echoes 
are observed beyond even that limit, as we now show. 

A theoretical decay curve may be calculated and com- 
pared to the experiments in Fig. 1, starting from a gen- 
eral spin Hamiltonian for ^^Si in doped Silicon. For ex- 
ample, for sample d, we have: 

Ti- = 'Hlab + Ti-^^Si-^^Si + 'W29 5j_3ip -I- 7^29 5j_e- , (1) 

where Tiiab includes the magnetic coupling of ^^Si spins to 
both the static laboratory field and the time-dependent 
tipping field produced by the rf pulses. Since ^^Si is fairly 
dilute (4.67% n.a.), 7^29 5^.29 5^ is just the direct dipolar 
coupling. The last two terms, Ti.29gi_3ip and Ti,29gi_f,- , 
play the role of the "bath" for the ^^Si spins, which pro- 
duce static magnetic shifts and determine Ti. In princi- 
ple, the dynamics of this bath might also affect our T2 
measurements. However, Fig. 1 shows that this is not 
the case, since samples a-d have nearly identical CPMG 
tails despite very different baths. This is strong empirical 
evidence that the ^^Si homonuclcar spin-spin coupling is 
sufRcient to describe the physics of all four samples (a-d), 
which greatly simplifies the model. Therefore, in the ro- 
tating frame, the secular part of Eq. (in the absence 
of rf pulses) is Hr , given by j9| : 




where il^ is the magnetic shift for spin i (relative to on- 
resonance spins), — - — — 3cos^6',;j] (^^7 is the 

gyromagnetic ratio for ^^Si), and bij = — y^- The vector 
between spins i and j , r^j , satisfies • z ~ cos 9ij . 

If some of the terms in Eq. ||2Jl are truncated, cor- 
responding to specific physical limits, then analytic so- 
lutions for the effect of various pulse sequences may be 
found using the product operator formalism III23. We 
start from the initial equilibrium density matrix: 

p(t = 0)cx J2 (3) 

i 

which assumes the conventional strong field and high 
temperature approximations 9]. 

For "unlike spins", where |ay| <ti Afi^ = jfi^ — 
we truncate the bij terms j9]. In this limit, the peak of 
the kth CPMG echo decays according to: 

{ly {k X TE)) = X: (0)f ff c^os (^M(h^\\ , (4) 
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FIG. 2: CPMG echo peaks at RT (triangles) and 4.2 K 
(squares) for sample d. The sets are scaled so that SEl agree. 
The solid line is the calculated decay from Fig. 1 (see text). 
While the qualitative temperature effect is clear, the non-ideal 
conditions of the 4.2 K data set (173y pulses and repetition 
time=100 sec ~ -^) prevent a quantitative assessment. 

which assumes the "infinite Hi limit" . Experimentally, 
^^^^sa22kHz, |^|<0.8 kHz, |^|<0.3 kHz for samples 
(a-c) and |§7l<3 kHz for sample (d). Eq. Q also de- 
scribes a free induction decay (FID) following a single 
90x pulse in another limit: all the bij terms are truncated 
and all fli = 0. Thus, the truncated dipolar decay of the 
CPMG echoes for the case of "unlike spins" is apparently 
unaffected by the 180y pulses, which flip all Iz^ — > {~Izi), 
leaving the sign of the bilinear terms unchanged. In 
order to compare Eq. to the data, we only need to 
have realistic values ofaij for our powder samples. To ob- 
tain these Qij, we built 20,000 "chunks" of the real Silicon 
lattice with arbitrary orientations, and determined the 
~80 nearest neighbors occupied according to the 4.67% 
n.a.. Averaging Eq. Q over all "chunks" yields the 
black curve shown in Fig. 1, which agrees remarkably 
well with the Hahn echo data points, but which fails to 
describe the measured CPMG echoes. 

When CPMG experiments are carried out in liquids, 
the well-known echo modulation due to J-coupling be- 
tween unlike spins can be effectively turned off |2ll |. 
if pulses are applied so frequently that i^^Jij and 
T^^AV,ij. Similarly, in our solid-state measurements, 
applying a CPMG sequence with frequent pulses (i.e., 
small TE) might push the system artificially into the "like 
spin" regime, where |ay|^Ar2y . In that limit, all the 
terms of Eq. Q should be retained. This precludes an 
analytic solution, but numerical calculations of (/y(i)) 
can be carried out for small numbers of spins, including 
the required ensemble averaging These calculations 
show that the initial decay of the CPMG echoes in that 
limit should be w | faster than Eq. I^J, which a gree s 
with the well-known second moment expressions |9|, |22| . 
Our data require another explanation. 

Empirically, the long tail induced by the CPMG se- 
quence has several interesting characteristics For 
example. Fig. 2 shows that the tail height, relative to 
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FIG. 3: CPMG echo trains at 4.2 K for sample d with TE of 
(a) 1.12 msec, (b) 2.65 msec, and (c) 11.23 msec. The solid 
line from Fig. 1 is scaled to intercept the first echo in each 
graph. The numbered echoes in (c-d) exhibit a pronounced 
even-odd asymmetry, which emerges for TE > T2jf^ . (d) 
shows the same effect in Si:Sb at RT with well-calibrated pulse 
angles, low repetition rates, and a narrow spectrum. 



the 1st echo, grows as the sample is cooled down to 4.2 
K. This result is only qualitative, since the tail height 
can also be changed by using a repetition time < Ti and 
by using tip angles slightly away from 180 degrees, which 
can be difficult to avoid at T=4.2 K. Still, even taking 
these factors into account, the tail appears to be more 
pronounced at low temperatures [l9| . 

To see if the long tail was due to some kind of multiple- 
pulse spin locking, we increased the interpulse spacing 
(TE), which led to another unexpected result. Fig- 
ure 3(a-c) shows data taken at 4.2 K for three differ- 
ent TE in sample d. The long tail persists even for 
TE> T2he- Interestingly, for large interpulse spacings, 
the odd-numbered echoes are much smaller than the 
even-numbered echoes. At RT, samples a-d exhibit the 
same even-odd asymmetry as TE is increased (e.g., as in 
Fig. 3(d)). 

This even-odd asymmetry leads to remarkable results 
as TE is increased still further. Figure 4 shows the FID 
and first two spin echoes acquired in a CPMG experiment 
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FIG. 4: The free induction decay (FID) and first two spin 
echoes (SEl, SE2) excited by a CPMG sequence at RT for 
sample d with TE of (a) 30 msec and (b) 60 msec. The insets 
show the narrow shape and the height of SE2 in comparison 
with SEl (the FID starts at 14600). At 60 msec, SE2 (a) is 
clearly different from SEl (b). Solid bars indicate pulses. 



with n = 2, for very long TE. In Fig. 4(a) 



TE 
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so that SEl is tiny relative to the FID. Surprisingly, SE2 
(at « 10.7) is nearly three times the height of SEl; 

SE2 is also narrower than SEl. In Fig. 4(b), TE is 
doubled, which pushes SEl into the noise, while SE2 is 
clearly visible, even though it occurs ^i.^xT2„g after the 
90x pulse. 

Since SE2 is the first echo to occur after three pulses, 
we decided to look for a contribution to the CPMG echoes 
that is reminiscent of a stimulated echo [lo| , using the se- 
quence: 90x-(^)-180F-TM-180y-DETECT, where TE 
and TM can be varied independently. Using this se- 
quence, we detect a conventional spin echo SE2 that 
peaks at total time 2 x TM, along with an "anomalous 
stimulated echo" (STE^) that peaks at TM-hTE. Figure 
5 shows the height of the STEa as either TM or TE is 
varied. There are several remarkable features of the data 
in Fig. 5: 1) we observe STE^, even for our best 180^ 
pulses, where there should be none, 2) they decay slowly 
as TE or TM are increased, 3) they appear to "start" at 
non-zero values at the left edge of Fig. 5, and 4) the data 
set has a larger scatter than expected, given the signal 
to noise of each individual data point. 

Given the results in Figs. 1-5, we have tried to min- 
imize the effects of nonidealities commonly reported in 
multiple pulse NMR (91,(13,113: (i) inhomogeneous Hi, 
(ii) finite-size JJi, (iii) a "spin locking" effect, and (iv) 
phase transients [ig. For (i), the results are unchanged 
if we use a tiny (~ 6% ) coil filling factor, or samples of 
very different skin depths. For (ii), the same effects are 
seen in all samples, even though i/i/FWHM changes by 
factor of 10. For (iii), similar results are obtained with 
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FIG. 5: "Anomalous stimulated echo" amplitudes at RT 
for sample d. Filled squares (TE~0.4 msec) are plotted 
vs. TE+TM. Empty circles (TM«10 msec) and triangles 
(TM«21 msec) are plotted vs TE. The solid line is from Fig. 
1. (Inset) The signal does not appear to grow from zero. 

an alternating phase Carr-Purcell sequence, where 180 
degree pulse phases alternate between"— X" and "X", 
even though the average Hi is quite different from that of 
CPMG. Finally, we expect that (iv) becomes less impor- 
tant as the number of pulses is reduced and their spacing 
is increased, so we don't see how this could explain the 
puzzling results of Figs. 4-5. 

Taken together, these results strongly suggest that the 
effects are due principally to the ^^Si homonuclear dipo- 
lar coupling. In that case, why is it so hard to find a 
quantitative explanation for the data? The form of Eq. 
(0) for a clean Silicon sample is one problem, since many 
spins may have |aij| ~ Afla, which make simulations 
[ijj particularly challenging [23- The dilution of the 
moments on the lattice could be another issue The 
strange features in Fig. 5 (in particular, point 3), and 
the narrowness of SE2 in Fig. 4 seem to be beyond the 
conventional theory of solid-state NMR. In recent NMR 
experiments j25j , large polarizations have produced mea- 
surable dipolar field effects, which led some to question 
the approximations underlying Eq. 0. While we don't 
have such large polarizations, the effects do appear to be 
more pronounced at low temperatures (Fig. 2). 

In the broader context of QC, the generic form of Eq. 
suggests that similar surprises may be found in other 
systems with small, long-range, qubit-qubit interactions, 
particularly when "bang-bang" control is used [2^. Un- 
derstanding these phenomena in Silicon may help to pre- 
vent similar surprises from imposing a performance limit 
on quantum computers in the future. 
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